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Introduction 


Prior Knowledge and Experience 


Using ArcGIS 9.3 or 10.0? 
— Workshop is for 10.0 


Prior Python experience? 


— I’m not assuming any 


Other programming experience? 


— I’m not assuming any 


Example 1 
* Script to copy all shapefiles in a folder into a geodatabase 


import arcpy 
from arcpy import env 
env.overwriteOutput = True 
env.workspace = "c:/workshop/ex01" 
fclist = arcpy.ListFeatureClasses () 
for fc in fclist: 
fcdesc = arcpy.Describe (fc) 
arcpy.CopyFeatures_management (fc, "c:/workshop/ex01/study.mdb/" 
+ fcdesc.basename) 


Example 2 


* Script tool to generate a k-nearest neighbor table 
e Runs an existing ArcGIS tool multiple times, writes the result 


import arcpy 

from arcpy import env 

env.overwriteoutput = True 

infc = arcpy.GetParameterAsText (0) 

output = arcpy.GetParameterAsText (1) 

k = arcpy.GetParameter (2) 

n= 1 

f = open (output, "w") 

while n <= k: 
result = arcpy.CalculateDistanceBand stats (infc, n) 
f.write(str(n) + " " + str(result[1])+ "\n") 
n=n+1 

f.close() 


Example 3 


* Script tool to run Huff model 
* Sophisticated analysis not available in ArcGIS 
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Example 3 


HuffModel. py - C:\Paul\Python\Examples\MarketAnalysisToolbox_update01262010\Script\Huf... E (©) 
Ele Edt Format Run Options Windows Help 
É HuffModel.py 
# Created: 4/13/2007 by Drew Flater 
Creating probability-based trade areas for retail stores 


# Import system modules 
import sys, string, arcgisscripting, os, traceback, shutil, re 


# Create the Geoprocessor object 
gp = arcgisscripting.create (93) 


# Set overwrite 
gp.overwriteoutput = 1 


def AddPrintMessage (msg, severity): 
print msg 
1f severity == 0: gp. AddMessage (msg) 
elif severity == 1: gp. AddWarning (msg) 
elif severity == 2: gp. AddError (msg) 


# Start traceback Try-Except statement: 
try: 
# Script parameters... 
stores = gp.getparameterastext (0) 
store_name = gp.getparameterastext (1) 
store attr = gp.getparameterastext (2) 
outfolder = gp.getparameterastext (3) 
fc name = gp.getparameterastext (4) 
studyarea = gp.getparameterastext (5) 
blockgroups = gp.getparameterastext (6) 





What is Python Scripting? 


* Add functionality to ArcGIS 
— Integrated into ArcGIS interface 
— Builds upon existing functionality 
— Automates repetitive tasks 
— Expands analysis options 


* Share new functionality 
— Script tools work just like regular tools 
— Can be integrated into models, tools 
— Easy to share with others (free) 


Why Python? 


Free, open source 

Object oriented 

Basic scripting AND complex object-oriented programming 
“Batteries included” 

Embraced by geospatial community, including ESRI 

Many libraries 


e. python 


&. python 


Python Community 


search 


Advanced Search 





ABOUT 
NEWS » 
DOCUMENTATION > 
DOWNLOAD » 
TÉ » 
COMMUNITY > 
FOUNDATION > 


CORE DEVELOPMENT > 


Help 
Package Index 
Quick Links (2.7.2) 

» Documentation 

» Windows Installer 
» Source Distribution 
Quick Links (3.2.2) 

» Documentation 

» Windows Installer 
» Source Distribution 
Python Jobs 

Python Merchandise 
Python Wiki 


Python Programming Language — Official Website 


Python is a programming language that lets you work more quickly and integrate 
your systems more effectively. You can learn to use Python and see almost 
immediate gains in productivity and lower maintenance costs. 


Python runs an Windows, I inuxitinix, Mac. OS X, and has been ported to the Java and 
NET virtual machines. 


Python is free to use, even for commercial products, because of its OSi-approved open 
source license: 


New to Python or choosing between Python 2 and Python 3? Read Python 2 or Python 3. 


The Python Software Foundation holds the intellectual propeny rights behind Python, 
underwrites the PyCon conference, and funds other projects in the Python community. 


Read more, -or- download Python now 

» PyArkansas 
The 4th annual PyArkansas will be held October 22, 2011 
Published: Wed, 21 September 2011, 11:37 +0200 


» PyGotham 
PyGotham will be held September 10-17, 2011. 


Published: Tue, 6 September 2011, 08:00 +0200 


http://www.python.org 


| wish there was Python 3 support in 


(enter PyPI package name) 


Tribon uses Python... 





joining users such as 
Rackspace. Industrial Light and 
Magic, AstraZeneca, Honeywell, 
‘and many others. 





hat they are sayin 
ITA Software: 





Python and ArcGIS 
e Python is the preferred scripting language for ArcGIS 


1. You can run Python from within ArcGIS 
— Python Window works like an interactive interpreter 

2. All tools in ArcToolbox can be accessed from Python 
— Import ArcPy to get full library of tools 

3. Python scripts can be made into tools 
— Extend functionality of ArcGIS 

4. Support for other scripting languages will go away 
— VBScript and JScript being replaced by Python 


Python Versions and ArcGIS 


Versions: 

— Current version of Python is 3.2.2 

— Python that works with ArcGIS 10.0 is 2.6.x 
Python that works with ArcGIS 10.1 is 2.7.x 
— Moveto Python 3.x likely only with ArcGIS 11 


ArcGIS only works with a specific version of Python: 
— Use the one that comes installed with ArcGIS 
— Don't install your own version of Python 


Installing Python 


Remove any existing installations of Python 


Install ArcGIS 10.0 
— Python 2.6.5 will be installed by default 


Install a Python editor 
Configure the editor to work with ArcGIS 


Note: You can run different versions of Python on one 
machine — however, a clean install of Python2.6.5 
with ArcGIS 10.0 is recommended 


Demo: Check ArcGIS and Python installation 


Fundamentals of Geoprocessing 
in ArcGIS 


Geoprocessing Tools 


Output Data 


( Input Data 


Lo a 





Tool Organization 


ArcToolbox 


681 & Analysis Tools 
& 


ata Interoperability Tools 
ata Management Tools 
& & Editing Tools 
E seocoding Tools 


inear Referencing Tools 
Multidimension Tools 


i & server Tools 
& & Spatial Analyst Tools 


Es Tracking Analyst Tools 





= 8) Analysis Tools 
2% Extract 


^. Table Select 
i$ & Overlay 
E & Proximity 
E & Statistics 


Tool Dialogs 


‘© Output Feature Class 





XY Tolerance (optional) 
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Tool Parameters 


* Parameters 


. Input Features 
— Required pem 
XY Tolerance (optional) 


— Optional 


e E r ro rs € Input Features 


C:\mydata\roads.shp 


A Output Feature Class 


* Warning Citat doo 


Environment Settings 


Environment Settings 


Y Workspace 

Y Output Coordinates 

Y Processing Extent 

Y XY Resolution and Tolerance 
Y M Values 

Y 2 Values 

Y Geodatabase 

Y Geodatabase Advanced 
Y Fields 

Y Random Numbers 

Y Cartography 

Y Coverage 

Y Raster Analysis 

Y Raster Storage 

Y Geostatistical Analysis 
Y Terrain Dataset 


Y TIN 





ok Cancel Show Help >> 





Geoprocessing Options 


Geoprocessing Options 


Overwrite the outputs of geoprocessing operations 
Log geoprocessing operations to a log file 


Background Processing 


Sto Notification J 
Appear for how long (seconds) 
Script Tool Editor/Debugger 


Editor: 





Debugger: 





Modetbuilder 


When connecting elements, display valid parameters when more than one is 
available. 


Results Management 
Keep results younger thar 


Display | Temporary Data 


[V] Add results of geoprocessing operations to the display 
[O results are temporary by default 





Demo: Geoprocessing Fundamentals 


Running Python Code 


Two ways to run Python Code 


1. Using an Interactive Interpreter 


— Code is executed directly line-by-line 


2. By running a script 
— Code saved in a .py file 


— Run from within a Python editor or directly from 
operating system 


Where to type and run Python code? 


1. Python window in ArcGIS 
— Built into any ArcGIS Desktop application 
— Good for testing code, very short scripts 


2. Python editor 
— IDLE installed by default 
— Many others, PythonWin is a good one to start 
— Good for more complex code, saving scripts 


Python Window in ArcGIS 


>>> print "Hello World" 
Hello World 


>>> 





Python Window in ArcGIS 


Works with current map document 
Interactive interpreter: 

— Executes code directly line-by-line 
Good for testing short code 

Code can be saved 

No error checking / debugging 


Python Editor - IDLE 


Python Shell 
Ele Edit shell Debug Options Windows Help 


Python 2.6.5 (r265:79096, Mar 19 2010, 21:48:26) [HSC v.1500 32 bit (Intel)] on win32 
Type "copyright", "credits" or "license(]" for more information. 


Personal firewall software may warn about the connection IDLE 
makes to its subprocess using this computer's internal loopback 
interface. This connection is not visible on any external 


interface and no data is sent to or received from the Internet. 


IDLE 2.6.5 
>>> 


Ln: 12 (col: 4 





Python Editor - PythonWin 


PythonWin 
File Edit view Tools Window Help 


OSU rr oo BARES 


Interactive Window 


PythonWin 2.6.5 (r265:79096, Mar 19 2010, 21:48:26) [MSC v.1500 32 bit (Intel)] on win32. 
Portions Copyright 1994-2008 Mark Hammond - see 'Help/About PythonWin' for further copyright information 
>>> Hello World 


£i hellowin. py 


print "Hello World" 


[script 'Cridatalhelowin,py' returned exit code O 





Python Editor 


Stand-alone — outside of ArcGIS 
Interactive interpreter: 

— Executes code directly line-by-line 

Save code as script files (.py) 

Good for organizing more complex code 


Demo: Running simple Python code 


Python Documentation 


Python Documentation 


& python = 





ABOUT Python Programming Language — Official Website 

NEWS : Python is a programming language that lets you work more quickly and integrate. 
DOCUMENTATION — > your systems more effectively. You can learn to use Python and see almost 

DOWNLOAD » immediate gains in productivity and lower maintenance costs. 1 wish there was Python 3 support in 
T z Python runs an Windows, ! intxitInbc, Mac OS X, and has been ported to the lava and 

COMMUNITY > NET virtual machines. (enter PyPI package name) 
FOUNDATION > 





Python is free to use, even for commercial products, because of its OSi-approved open 
source license: 


CORE DEVELOPMENT > 








Help Tribon uses Python... 
Package Index New to Python or choosing between Python 2 and Python 3? Read Python 2 or Python 3. 
Quick Links (27.2) 
+ Dobumentation ‘The Python Software Foundation holds the intellectual property rights behind Python, 
» Windows Installer underwrites the PyCon conference, and funds other projects in the Python community. 
> Source Distribution 
itu na Read more, -or- download Python now Joining users such as 
RE oe Rackspace. Industrial Light and 
E Cani Magie, AstraZeneca, Honeywell 
» Windows Installer The 4th annual PyArkansas will be held October 22, 2011 
he and many others. 
as Published: Wed, 21 September 2011, 11:37 +0200 
Python Jobs +PyGatham 
Python Merchandise FyGotham will be held September 16-17, 2011 hat they are sayin 
Python wish Published. Tue, 6 September 2011, 08.00 +0200 MELIA 


http://www.python.org 


Python Documentation 


Version specific! 


Python modules | index 
Download Python v2.6.7 documentation 
d E Welcomel This is the documentation for Python 2.6.7, last updated Jun 03, 2011 


Parts of the documentation: 


What's new in Python 2.6? Extending and Embedding 


or all "What's new" documents since 2.0 tutorial for C/C++ programmers 
Tutorial Python/C API 
start here reference for C/C++ programmers 





Using Python 


how to use Python on different platforms 


Library Reference 
keep this under your pillow 


Language Reference 
describes syntax and language elements 


Python HOWTOs 
in-depth documents on specific topics 


http://docs.python.org 


Installing Python Modules 


information for installers & sys-admins 


Distributing Python Modules 
sharing modules with others 


Documenting Python 
guide for documentation authors 


FAQs 


frequently asked questions (with answers!) 


Python Beginners Guide 


® python 


FRONTPAGE 
RECENTCHANGES 
FINDPAGE 
HELPCONTENTS 
BeGInuersGUIDE 
Page 

» Immutable Page 
» Info 

» Attachments 

» More Actions: 
User 


» Login 





Search titles | text 


» BeginnersGuide 





» BeginnersGuide 


Beginner's Guide to Python 


New to programming? Python is free, and easy to learn if you know where to start! This guide will help you to get 
started quickly. 


Chinese Translation 
New to Python? 
Read BeginnersGuide/Overview for a short explanation of what Python is. 
Getting Python 


Next, install the Python interpreter on your computer. This is the program that reads Python programs and carries 
out their instructions; you need it before you can do any Python programming. 


There are currently two major versions of Python available: Python 2 and Python 3. The Python2orPython3 page 
provides advice on how to decide which one will best suit your needs. At the time of writing (21 Jun 2010), the rest 
of this page assumes you've decided to use Python 2. 


See BeginnersGuide/Download for instructions for downloading the correct version of Python. 


At some stage, you'll want to edit and save your program code. Take a look at HowToEditPythonCode for some 
advice and recommendations. 


http://wiki.python.org/moin/BeginnersGuide 


Python Books 


From Novice to Professional 





None of these books including anything on ArcGIS or geoprocessing! 


Python Language Fundamentals 


Python Data Types 


Number (integer and float) 
String 

List 

Tuple 

Dictionary 


Strings, lists and tuples are sequences 
Strings, numbers and tuples are immutable 
List and dictionaries are mutable 


Numbers 


* Integers 
— Whole number, i.e. no decimals 
— e.g. -34 
e Floats 
— Decimal point 
— e.g. -34.8307 


Numerical Operators 





Operator Description Integer Floating-point 
Example Result Example Result 
Multiplication 9*2 18 9*20 18.0 
/ Division 9/2 4 9/2.0 4.5 
% Modulus 9% 2 1 9% 2.0 1.0 
+ Addition 9+2 11 9+2.0 11.0 


- Subtraction 9-2 7 9-2.0 7.0 





Demo: Numerical Operators 


Strings 


* Aset of characters surrounded by quotes is called a 
string literal 


e To create a string variable, assign a string literal to it 


>>> mytext = "Crime hotspot maps are cool." 
>>> print mytext 


Crime hotspot maps are cool. 


Quotes in Python 


* In Python single and double quotes are the same 
* "NIJ" isthesameas 'NIJ' 


>>> print "I said: 'Let's go!'" 


* Quotes in Python are straight-up 
* "text" or 'text', not "text" or ‘text’ 


* Be aware of copy/paste and auto-formatting 


Variables 


* Python scripts use variables to store information 
* To work with variables use an assignment statement 


>>> x = 17 
>>> x * 2 
34 


Variables 


* Python uses dynamic assignment 


>>> x = 17 
>>> type (x) 
«type 'int'» 
>>> x = "GIS" 
>>> type (x) 
<type 'str'> 


* No need to declare variables 
* Value defines the type 


Variable Names 


* Rules 
— Letters, digits and underscores 
— Cannot start with a digit 
— Don't use keywords (print, import, etc.) 


* Recommendations 
— Be descriptive (count instead of c) 
— Keep it short (count instead of count of records) 


— Follow convention: all lowercase, use underscores 


Statement and Expressions 


* A Python expression is a value 


>>> 2 * 17 
34 


* A Python statement is an instruction to do something 


>>> x = 2 * 17 


Working with Strings 


e Concatenate strings 


>>> x = RGS 
>>> y = "I" 
>>> 2 = WOW 


>>> print x + y +z 
GIS 


Converting to String 


>>> temp = 100 
>>> print "The temperature is " + temp + " degrees" 


TypeError: cannot concatenate 'str' and 'int' objects 


>>print "The temperature is " + str(temp) + " degrees" 


* Converting the value of a variable from one type to another is 
known as casting 


Lists 


* A Python list is an ordered set of items 


* The list of items is surrounded by square brackets [ ], and the 
items are separated by commas (,) 


* Items can consist of numbers, strings and other data types 


mylist = [1, 2, 4, 8, 16, 32] 
mywords ["jpg", "bmp", "tif"] 


* Lists are very widely used in geoprocessing: 
— e.g. list of feature classes, list of records, list of fields, etc. 


Python Functions 


* A function carries out a certain action 
* Python has many built-in functions 


«function» («arguments») 


>> pow(2,3) 
8 


* Using a function is referred to as calling a function 
* Additional functions can be accessed using modules 


Python Methods 


* Amethod is a function that is closely coupled to some object 
<object>.<method> («arguments») 
>>> topic = "Crime Mapping" 


>>> topic.count ("i") 
2 


* Many of Python’s data types have methods 


String Indexing 


* Python strings have an index positioning system 


>>> mystring = "Crime Mapping" 
>>> mystring[0] 

Ecos 

>>> mystring[-1] 

tg! 


* Strings can be sliced into smaller strings using slicing 


>>> mystring[0:5] 
'Crime' 


Working with List 


+ Python lists have an index positioning system 


>>> crimes = ["arson", "burglary", "robbery"] 
>>> cities[1] 


'burglary' 
* There are many list methods 


>>> crimes.append("homicide") 
>>> crimes.remove ("arson") 
>>> crimes 


['burglary', 'robbery', 'homicide'] 


Working with Pathnames 


* Pathnames are critical when writing scripts: 
— Example workspace: c:\data\results 
— Example shapefile: c:\data\results\streams.shp 
* In Python a backslash (1) is an escape character 


e Pathnames in Python should therefore look like one 
of the following 


"c:/data" (mm 


"c:\\data" 
r"c:\data" (raw string) 


Python Modules 


e Modules are like extensions that can be imported into Python 
to extend its capabilities 


>>> import time 


e A typical module contains a number of specialized functions 
which can be called once the module has been imported 


<module>.<function> 


>>> time.localtime() 


Conditional Statements 


* Branching can be used to control workflow 


import random 

x = random.randint (0, 6) 
print x 

if x == 6: m 


print = "You win!“ 


* Syntax: keyword i £, followed by a condition, followed by (:) 


Indentation in Python 


Indented code is referred to as a block 
Use tabs or spaces — be consistent 
Recommended: 4 spaces 


Tip: be careful with copy/paste from other 
applications 


More Conditional Statements 


* Useofelif andelse is optional 


import random 
x = random.randint (0,6) 
print x 
if x == 

print "You win!" 
elif == 5: 

print "Try again!" 
else: 

print "You lose!" 


Loop Structures: While 


* Loop structures allow you to repeat a certain part of your 
code 


e Awhile loop repeats until a particular condition is reached 


i=0 

while i <= 10: 
print i 
i += 1 


* The while statement uses a sentry variable in the exit 
condition 


Loop Structures: For 


* A for loop repeats a block of code for each element of a 
sequence 


mylist = ["A", "Bm, nou "D"] 
for letter in mylist: 


print letter 


* Inthe example, letter is the name of a variable and for each 
iteration of the loop this varaible is assigned a different value 


ArcPy: Geoprocessing using Python 


What is ArcPy? 


* ArcPy was introduced with ArcGIS 10.0 


* ArcPy is a collection of modules, classes and 
functions which give access to all the geoprocessing 
tools in ArcGIS from within Python 


* Most geoprocessing scripts will start with: 
import arcpy 


* Note: ArcPy replaces the older arcgisscripting 
module 


Setting Current Workspace 


* After importing ArcPy, most scripts start with setting a 
workspace to retrieve and store files 


import arcpy 


arcpy.env.workspace = "c:/workshop" 


* Inthe code above env is a class and workspace isa 
property of this class 


arcpy.<class>.<property> 


Using Tools 


* ArcPy gives you access to all tools in ArcToolbox 
* All tools are provided as functions 


arcpy.«toolname toolboxalias» («parameters») 


* Example: 


import arcpy 
arcpy.env.workspace - "c:/data" 


arcpy.Clip analysis("streams.shp", "study.shp", "result.shp") 


Tool Parameters 


A good understanding of tool parameters is essential 
Parameters have properties: 

— Name 

— Type (feature class, integer, etc.) 

— Direction (input or output) 

— Required or optional 


Example: Clip tool 


Parameter Explanation Data Type 

in features The features to be clipped. Feature Layer 
cip Joras The features used to clip the input features. Feet Layer 
out feature dass The feature dass to bescrested: Feature Class 


cluster. tolerance 


The minimum distance separating all feature coordinates (nodes and Linear unit 
(Optional) vertices) as well as the distance a coordinate can move in X or Y (or 
both). Set the value to be higher for data with less coordinate 


accuracy and lower for data with extremely high accuracy. 


Tool Syntax 


Tool dialog: = 


Clip Features 











‘Output Feature Class 








XY Tolerance (optional 














(x — Cea] [Erin] [E] 








Python syntax: Clip analysis(in features, clip features, 
out feature class, 
(cluster tolerance]) 


Example: Clip analysis ("streams.shp","study.shp", 
"result.shp") 


Optional Parameters 


* Required tool parameters are listed first 


* Optional tool parameters can be left out 
— But what if some need to be set? 


Buffer_analysis (in_features, out_feature_class 
buffer distance or field, (line side), (line end type), 
(dissolve option), (dissolve field)) 


arcpy.Buffer analysis("roads", "buffer", "100 METERS", "", 
"", "LIST", "Code") 


arcpy.Buffer analysis("roads", "buffer", "100 METERS", 
dissolve option-LIST, dissolve field-Code) 


Hard-coded Parameters 


* Consider the example 


import arcpy 
arcpy.env.workspace = "c:/data" 


arcpy.Clip_analysis("streams.shp", "study.shp", "result.shp") 


* How can we make this code more usable? 


Using Variables for Parameters 


import arcpy 


arcpy.env.workspace = "c:/data" 
infc = "streams.shp" 
clipfc = "study.shp" 
outfc = "result.shp" 


arcpy.Clip analysis (infc, clipfc, outfc) 


Variables Provided by a User 


import arcpy 

infc = arcpy.GetParameterAsText (0) 
clipfc = arcpy.GetParameterAsText (1) 
outfc = arcpy.GetParameterAsText (2) 


arcpy.Clip analysis (infc, clipfc, outfc) 


Result Objects 


* ArcPy returns the output of a tool as a Result object 


import arcpy 
arcpy.env.workspace = "c:/data" 
myresult = arcpy.Clip analysis("streams.shp","study.shp","result.shp") 


print myresult 


* This will print the path to the output dataset 


c:/data/result.shp 


Multiple Operations using Result Objects 


* Result objects can be used as the input into another 
function 


import arcpy 

arcpy.env.workspace = "c:/data/study.gdb" 

buffer = arcpy.Buffer_analysis("str","str_buf","100 METERS" 
count = arcpy.GetCount_management (buffer) 


print count 


* This allows complex geoprocessing operations 


ArcPy Classes 


* Some tool parameters are complicated/detailed 
— e.g. coordinate system 


* ArcPy classes are used to work with these parameters 
— Classes are used to create objects 
— Classes have properties and methods 


* General syntax 


arcpy.<classname> (<parameters>) 


ArcPy Classes: Example 


* The following is an example of the contents of a .prj file 


B. streams. prj - Notepad 
File Edit Format View Help 


|PRO3CS["NAD_1983_statePlane_Florida_East_FIPS_0901_Feet" ,GEOGCS [ GCS 
North american 1983", DATUM["D North. American 1983", SPHEROID["GRS_19 
80", 6378137. 0,298.257222101]], PRIMEM[ “Greenwich”, O. òl; UNIT "Degree", 





9. 0174532025199433]], PROJECTION” Transverse Mercator’ ER PARAMETER ["Fa 
a * PARAMETER False.) Northing", 0.0], PARA 


PARAMETER [Scale Factor", 0. 999941176 
4705882], PARAMETER ("| Latitude oF. origin", 24. 33333333333333], UNIT[" Foo 
tus" *, 0, 3048006096012192]] 





* To avoid having to work with this actual string, we can use a 
SpatialReference class 


ArcPy Classes: Example 


* The following example creates a spatial reference object 
based on an existing .prj file - properties of this object can 
then be used 


import arcpy 

prjfile = "c:/data/streams.prj" 

spatialref = arcpy.SpatialReference (prjfile) 
myref = spatialRef.name 

print myRef 


* This will print 


NAD 1983 StatePlane Florida East FIPS 0901 Feet 





ArcPy Classes: Example 


* The following example creates a spatial reference 
object and use this to define the coordinate system 
of a new feature class 


import arcpy 


out_path = "c:/data" 
out_name = "lines.shp" 
prjfile = "c:/data/streams.prj" 


spatialref = arcpy.SpatialReference (prjfile) 
arcpy.CreateFeatureclass management (out path, out name, 
"POLYLINE", "", nn, vn, spatialref) 


ArcPy Functions 


* All geoprocessing tools are ArcPy functions 
* Additional ArcPy functions: 

— listing data 

— Retrieving and setting properties 

— Many more... 


* General syntax 


arcpy.«functionname» («arguments») 


ArcPy Functions 


Cursors 

Describing data 

Environment and settings 
Fields 

General 

General data functions 
Getting and setting parameters 
Licensing and installation 
Listing data 

Messaging and error handling 
Progress dialog 

Tools and toolboxes 


Describing and Listing Data 


Describing Data 


* The Describe function is used to determine properties of dataset 
* General syntax 


import arcpy 
<variable> = arcpy.Describe (<input dataset>) 


* Example: 
import arcpy 


desc = arcpy.Describe("c:/data/streams.shp") 
print desc.shapeType 


Describing Data: Example 


import arcpy 


arcpy.env.workspace = "c:/data" 
infc = "streams.shp" 

clipfc = "study.shp" 

outfc = "streams clip.shp" 


desc = arcpy.Describe (clipfc) 
type = desc.shapeType 
if type == "Polygon": 
arcpy.Clip analysis (infc, clipfc, outfc) 
else: 
print "The clip features are not polygons." 


Listing Data 


* Listing data is very common 


* Several different list functions in ArcPy 
— ListFields 
— Listlndexes 
— ListDataset 
— ListFeatureClasses 
— ListFiles 
— ListRasters 
— ListTables 
— ListWorkspaces 
ListVersions 


. Similar logic: 


— Create a list 
— Iterate over the list using a for loop 


Listing Feature Classes 


* TheListFeatureClasses function returns a list of 
feature classes in the current workspace 


* General syntax: 


ListFeatureClasses ((wild card), (feature type), 
(feature dataset)) 


* Example: 


import arcpy 

from arcpy import env 
env.workspace - "c:/data" 

fclist = arcpy.ListFeatureClasses() 


Listing Feature Classes 
* No filtering: 
fclist = arcpy.ListFeatureClasses () 
* Filtering based on wild card 
fclist = arcpy.ListFeatureClasses ("w*") 
* Filtering based on feature type 


fclist = arcpy.ListFeatureClasses("", "point") 


Listing Fields 


* The ListFields function lists the fields in a feature class or 
table in a specified dataset. 


* General syntax: 

ListFields (dataset, (wild card), (field _type)) 
* Example 

import arcpy 


arcpy.env.workspace - "c:/data" 
fieldlist = arcpy.ListFields ("roads.shp") 


Using Lists in for loops 


* The following script creates a list of fields of type 
String and determines for each text field what the 
length of the field is 


import arcpy 
arcpy.env.workspace = "c:/data" 
fieldlist = arcpy.ListFields("roads.shp", "", 
"String") 
for field in fieldlist: 
print field.name + " " + str(field.length) 


Using Lists in for loops 


* The following script creates a list of TIFF files and 
iterates through each file in the list to build 
pyramids 


import arcpy 

from arcpy import env 

env.workspace = "c:/data" 

tifflist = arcpy.ListRasters("","TIF") 

for tiff in tifflist: 
arcpy.BuildParamids management (tiff) 


Creating Custom Tools 


Ways to Execute a Script 


1. Asa stand-alone script 
— The script is executed from the operating system or from 
within a Python editor such as PythonWin 
— When using ArcPy, ArcGIS needs to be installed and 
licensed 
— No ArcGIS Desktop application needs to be open 


2. Asa script tool within ArcGIS 

Atool dialog is created to execute the script 
Script tool looks like any other tool in ArcToolbox 
Tool execution is controlled from ArcGIS Desktop 
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Why Create Script Tools? 


Tool dialog makes it easier to use 

Tool dialog validates user inputs 
Becomes part of all geoprocessing 
Environment settings are passed on 
Writes messages to the Results window 
Easier to share 

Does not require user to know Python 
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Steps to Create Script Tools 


Create a Python script (.py) 

Create a custom Toolbox (.tbx) 

Add a tool to the Toolbox using Add Script 
Modify the script with inputs and outputs 


Example Script: Hardcoded Variables 


import arcpy 
from arcpy import env 
env.overwriteoutput = True 





f = open (output, "w") 

while n <= k: 
result = arcpy.CalculateDistanceBand stats (infc, n) 
f.write(str(n) + " " + str(result[1])+ "\n") 
n=n+ 1 

f.close() 
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Example Script: User Provided Parameters 


import arcpy 
from arcpy import env 
env.overwriteoutput = True 


f = open (output, "w") 
while n <= k: 


result = arcpy.CalculateDistanceBand stats (infc, n) 
f.write(str(n) +" "+ str(result[1])+ "\n") 
n=n+ 1 


f.close() 


More ArcPy Functionality 


More ArcPy Functionality 


Cursors to work with rows and geometry 

— Retrieve, edit, create 

arcpy . sa module to work with rasters 
arcpy .mapping module for map automation 
Creating custom functions and classes 


Resources for Python Scripting in ArcGIS 
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ArcPy is a site-package that builds on (and is a successor to) the successful arcgisscripting 
module. Its goal is to create the cornerstone for a useful and productive way to perform 
geographic data analysis, data conversion, data management, and map automation with 
Python. 

This package provides a rich and native Python experience offering code completion (type a 
keyword and a dot ta get a pop-up list of properties and methods supported by that 
keyword; select one to insert it) and reference documentation tor each function, module, 
and class. 

The additional power of using ArcPy within Python is the fact that Python is 3 general 
purpose programming language. tt is interpreted and dynamically typed and is suited for 
interactive work and quick prototyping of one-off programs known as scripts while being 
powerful enough to write large applications in. ArcGIS applications written with ArcPy 
benefit from the develanment af additional modules in numerous niches af Python hy GIS 
professionals and programmers from many different disciplines. 


General Help 
python provides the facility of doiusêntatisa ‘ating: The functions and classes availble 
in ArcPy use this method for the package documentation. One method for reading these 


messages and getting help is by using the command ne1p provided by Python. Running 
the command with an argument displays the calling signature and the documentation 


string of the object. 


>>> import arcpy 
>>> help(arcpy) 
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Description 


AL ArcGIS Desktop 10, Python scripting is tightly integrated into ArcMap and ArcCatalog, allowing you to create 
and automate GIS workflows quickly and easily. This course introduces Python scripting in ArcGIS Desktop and 
shows how you can use scripts to increase productivity and the quality of ycur maps and data. The 
presentation covers how to use the new ArcPy mapping module to manipulate map documents and layers. 


Who Should Attend 


GIS analysts, specialists, and other experienced ArcGIS users who want to automate complex tasks and 
common nroredures. 


Goals 
After completing this course, you wil be able to 


+ Create basic Python scripts using correct syntax. 
+ Write and run scripts in ArcMap using the Python window. 

+ Use Python in the Field Calculator. 

+ Create script tools to automate geoprocessing operations. 
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This is a Free Course 


Questions? 
Contact us via e-mail or call 
toll-free at 888-377-4575, 
select option 3, between 
3:00 AM and 5:00 PM (Pacific 
Time). 
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Spatial analysis is one of the more interesting anc 
remarkable aspects of GIS. Using spatial analysis, you can 
combine information from mary independent sources and 
derive a new set of information (rasults)—by applying a 
large, rich, and sophisticated set of spatial operators. 
GIA nmfassinnals use gennrneessing rn pragram their mun 
* ideas in order to derive thasa analytical rasutts, In tum, 
these results are applied to a wide variety of problems, For 
examole, hera, geoprocessing is used to identify suitable sites for parks. The 
result is a dataset of potential park sites for further evaluation, Site selection 
logic is used to find areas that ere close to where people Ive but are not too 
close to existing parks. 
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Beyond ArcGIS 


Using PySAL for Spatial Analysis 
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PySAL 


* Python library of spatial analysis methods 
* ESDA, spatial statistics, geostatistics 
* Growing and expandable 


Using R for Spatial Analysis 


Open source language for data analysis 

Libraries have been developed for spatial methods 
Large and active user community 

Growing and expandable 


ArcGIS and R 
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Evaluating R Statements 
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Print CHelo morid"") 
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Speed <- cars$speed 
Distance <- carssdist 
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Concluding Remarks 


Python is a relatively easy to learn language 
ArcGIS is becoming more “Pythonesque” 


Creating time-savings scripts for repetitive tasks does 
not take a lot of code 


Easy to share script tools 
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